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Abstract 

We describe a universality class for the transitions of a generalized Polya urn by studying the 
asymptotic behavior of the normalized correlation function C[t) using finite-size scaling analy¬ 
sis. X(1),X(2), ■ • • are the successive additions of a red (blue) ball ( X(t ) = 1(0)) at stage t 
and C(t) = Cov(X(l),X(t + l))/Var(X(l)). Furthermore, z(t) = X(s)/t represents the 

successive proportions of red balls in an urn to which, at the t + 1 -th stage, a red ball is added 
( X(t + 1) = 1) with probability q(z{t)) = (tanh[J(2^(t) — 1) + h] + l)/2, J > 0, and a blue ball is 
added (X(t + 1) = 0) with probability 1 — q(z(t)). A boundary (J c (h), h ) exists in the (J, h ) plane 
between a region with one stable fixed point and another region with two stable fixed points for 
q(z). C(t) ~ c + c' ■ t l ~ l with c = 0 (> 0) for J < J c (J > J c ), and l is the (larger) value of the 
slope(s) of q(z) at the stable fixed point(s). On the boundary J = J c (h), C(t ) ~ c + c' ■ (lnt) -a, 
and c = 0 (c > 0), a! = 1/2 (1) for h = 0 (h 7 ^ 0). The system shows a continuous phase transition 
for h = 0 and C(t ) behaves as C{t) ~ (lnt) -o/ g((l — l)\nt) with a universal function g(x) and a 
length scale 1/(1 — l ) with respect to Inf. [3 = v\\ ■ a' holds with (3 = 1/2 and v\\ = 1. 

PACS numbers: 05.70.Fh,89.65.Gh 
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I. INTRODUCTION 


The contagion process is one of the most studied topics in statistical physics and has 
attracted the attention of many researchers from various disciplines jl]-8] . Polya urn is one 


of the simplest models for this process [9-llJ. In this model, an urn consists of t balls, where 
the proportion of red balls is z(t ) G (0,1) and the rest of the balls are blue. The probability 
of a new red ball being added to the urn is z(t ), while it is 1 — z(t ) for a new blue ball; the 
proportion of red balls then becomes z(t + 1). This procedure is iterative, which produces a 
sequence of proportions z(t 0 ), z(t 0 + 1), z(t 0 + 2), • • •, where the urn contained t 0 ■ z(t 0 ) red 
balls at t = to- The limit value lim^oo z(t ) obeys a beta distribution with shape parameters 
a = t 0 ■ z(t 0 ) and j3 = to ■ z(t 0 ). 

As the Polya urn process is very simple and there are many reinforcement phenomena in 
nature and the social environment, many variants of the process have been proposed, referred 
to as generalized Polya urn processes[12]. In the nonlinear generalizations of this model, a 
continuous function q : [0,1] — > [0,1] determines the probability q(z(t)) of a red ball being 


ded at stage t + 1. This nonlinear version is referred to as a nonlinear Polya process 
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ad 

Il4] . In contrast to the original linear model, the nonlinear model can have many isolated 
stable states. The fixed point z * of q(z), where q(z*) = z *, is (un)stable if (z — z*)(q(z) — z) 
is negative (positive) for all z in the vicinity of z* Q. «. , referred to a S d_irr g 
(upcrossing) as the graph y = f(z) crosses the curve y — z in the downwards (upwards) 
direction. The slope of q(z) at z* is smaller (larger) than one when z* is downcrossing 
(upcrossing). When q(z) touches the diagonal q — z in the (z,q) plane at z t , a point that 


is referred to as the touchpoint, the stability of z t depends on t 
slope of q(z) and the diagonal z in the left neighborhood of z t 


re difference between the 
. If it is less (more) than 


1/2, Zt is (un)stable. The multiplicity of the stable states provides a convenient picture 


that explains the lock-in phenomena in the technology and product adoption processes 


id- 


Suppose two selectively neutral technologies enter the market at the same time. Because 
economies of scale play the role of an externality that persuade a new consumer to buy the 
dominant technology, determining which technology to buy depends on the proportion of 
each technology possessed by previous consumers. If the dependence is described by the 
non-linear function q(z), the technology adoption process is described by a non-linear Polya 
urn. An S-shaped q(z) function with two stable fixed points suggests the random monopoly 
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formed when one technology dominating over the other depends on chance fluctuations at 
the start [? ]. 

Information cascade provides a good experimental setup for verifying theoretical predic¬ 
tions 
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19]. Here, participants sequentially answer questions with two possible choices 


(two-choice questions). In addition to their own information or knowledge, they refer to 
social information about the number of previous subjects that chose each option. In an 


experiment where subjects answer general knowledge two-choice questions, it is possib 


change the number of stable states by controlling the difficulty of the question |18l. [19]. A 


e to 


subject that knows the answer to a question chooses the correct answer with a probability 
of 1. A subject who does not know the answer tends to choose the majority choice. By 
changing the difficulty of the question, an experimenter can control the ratio p of the latter, 
no-knowledge, subject. We denote the probability that the no-knowledge subject choose the 
correct answer as qh{z), when the ratio of correct choices among previous subjects is z. The 
probability that a subject choose the correct answer is then q(z) = (1 — p) ■ 1 + p ■ qh(z). 
The sequential voting process in the experiment is described by a non-linear Polya urn. It 
was shown that q(z) has one (two) stable fixed point(s) for p < p c (> p c ). 

If there is only one stable fixed point, z(t) converges to it. In the case of multiple stable 
fixed points, the stable fixed point to which z converges will be random]13], By controlling 
the model parameters, the number of stable fixed points can be changed, which induces a 
non-equilibrium transition. In an exactly solvable case, where q(z) is a combination of the 
constant g* G (1/2,1] and the Heaviside function 9(z— 1/2), i.e., q(z) = (1— p)-q if +p-6(z— 1/2) 
with a correlation control parameter p G [0,1], q(z) touches the diagonal at Zt — 1/2 for 
p = p c = 1 — 1/2g*. For p < p c , there is a unique stable state at z+ — (1 — p)q * + p. For 
p > p c , there are two stable states at z± = (1 — p)q ±p. Because the slope of q(z) in the left 
neighborhood of z t is 0, the touchpoint at z t is unstable, and z(t) converges to the stable 
fixed point at z + for p = p c [20]. 

The probability of convergence to a stable fixed point depends strongly on the color of 
the first ball when there are multiple stable states. If the color is red and z( 1) = 1 (blue 
and z(l) = 0), the probability of convergence to a larger stable fixed point becomes higher 
(lower). The difference in the probabilities is given by the limit value c of the normalized 
correlation function C(t) between the color of the first ball and that of the t + 1-th ball. 
Furthermore, c plays the role of the order parameter for the phase transition. In the afore- 
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mentioned exactly solvable model, which we refer to as the ’’digital” model, C(t) at p = p c 
shows a power law dependence on t as C(t) oc t~ a with a = 1/2. The order parameter be¬ 
haves as c oc (p — PcY with f3 = 1. In addition, C{t) obeys the scaling form C(t ) oc t~ a g(t/^) 
near p c with a universal function g and correlation length £. £ diverges as £ oc \p — p c I - " 11 
with v\\ — 2 2 ill. The scaling relation /? = v\\ ■ a holds as in the absorbing states phase 


transition 


22 


23]. 


In this work, we use finite-size scaling (FSS) analysis in order to study the asymptotic 
behavior of the correlation function for generalized Polya urns. We adopt a logistic-type 
model q(z) = (tanh[J(2z — 1) + h\ + l)/2 with two parameters J and h. Here, J is the 
parameter of the strength of the correlation, and h is the parameter of the as ym metry. The 
motivation for adopting this model was derived from experimental findings 18), [l9j]. With 


this choice, there is a threshold value J c {h), and at J = J c (h), q(z) becomes tangential to 
the diagonal at z t . From the above discussion, the touchpoint at z t is stable, which differs 
from the digital model. There are two stable states at z t and z + for /i / 0. If the order 
parameter c takes a positive value at J = J c (h ), the phase transition becomes discontinuous. 
For h — 0, the touchpoint at z t = 1/2 is the unique stable state, c should be equal to zero 
and the phase transition becomes continuous. We also clarify the universality class of the 
continuous transition. 


The remainder of the paper is organized as follows. Section [TT] introduces the model. In 
section nm we discuss the asymptotic behavior of C(t) for J ^ J c (h ) using previous results, 
and propose the asymptotic form C(t) — c-t-c^lnt)^ for J = J c (h). In section HVl we study 
the FSS relations of the system. Section|V]is devoted to the FSS study of C(t) for J = J c (h). 
We show that c > 0 (c = 0) and a' = 1 (a' = 1/2) for h > 0 [h — 0). The system shows 
a continuous phase transition for h = 0. In section IVII we study the universality class of 
the continuous transition. We show that C(t) oc (lnt) _Q/ g(lnt/£) with a universal function 
g{x) and a length scale £. We define the critical exponents /3 and v\\ as c oc (J — J C )P and 
£ oc | J — J c \~ Uu , respectively. Using the scaling relation /3 = v\\ ■ a' with v\\ = 1, we obtain 
/3 = 1/2. Section IVIII provides our summary and further comments. In Appendix A, we 
derive the explicit form of g(x) for the digital model. 
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II. MODEL 


We define the stochastic process X(t) G {0, l},f G {1,2, •• • , T}, where the probability 
that X(t) takes value of 1 is given by a function q(z) of the proportion z(t — 1) of the 
variables X(l), • • • ,X(t — 1) that are equal to 1. 


q(z) = Pr(X(f) = l\z(t — 1) — z) — - (tanh [J( 2 • z — 1) + h] + 1) , 

1 t i 

= t E X(s) for t > 0, and z(0) = -. (1) 


The choice of q(z) is arbitrary, and we adopt the above form, which is familiar in the held 
of physics 24,125]. The fixed point of q(z) is a solution to q(z) = z. Using the mapping 
m = 2z — 1, we obtain the self-consistent equation m = tanh( J ■ m + h) for the magnetization 
m in the mean-held Ising model. Here, we consider only the case for which J > 0. Because 
of the symmetry under (.X , h) GG (1 — X , —h), we also assume that h > 0. 



FIG. 1. Plot of J c (h) in (J,h). 

The number of hxed points for q(z) depends on ( J,h ). There is a threshold value J = 
J c (h) as a function of h (Fig. [TJ). For J < J c (h), there is only one hxed point at z = z + 
(FigJ2{a)). With increasing J, q(z) becomes tangential to the diagonal z at z t for J = J c (h). 
For h > 0, z t 7 ^ z + , and both z t and z + are stable (FigJ3(b)). For h — 0, z t — z + , and 
it is also stable (Figj3](a)). In both cases, the slope of the curve at z t is equal to one. For 
J > J c (h), there are three hxed points, and we denote them as < z u < z + ; z± is stable, 
and z u is unstable (Figj2{b)). We denote the value q(z±), q(z t ) as q±, q t and the slope of q(z) 
at z±,Zt as l± = q'(z±),l t — q'( z t ) — 1- As z± is stable and downcrossing, l± < 1 - 
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FIG. 2. Plot of q(z) vs. z for h > 0. The intersection between y = q(z) and y = z is the fixed 
point of q(z). (a) J < J c {h), with one fixed point at z = z + . (b) J > J c (h), with three fixed points 
at z £ {z-,z u ,z + }. 




FIG. 3. Plot of q{z) vs. z for J = J c {h) (thin solid line) and J > J c (h) (thick solid line) for (a) 
h = 0 and (b) h > 0. 


We note a crucial difference between h — 0 and h > 0. For h — 0, the touchpoint at 
z t — 1/2 for ,J = J c (0) = 1 coincides with the stable fixed point at z + — 1/2 for J < J c (0). 
It splits into the two stable hxed points at z — z± for J > J c (0). (Figj3](a)). z± continuously 
moves away from z t and z + — Z- oc (J — Jc) 1 ^ 2 for | J — J c \ « 1 as in the case of the mean 
field Ising model. On the other hand, for h > 0, the touchpoint z t appears at a different 
position from z + for J = J c {h ) (Figj3])b)). As J increases from J c (h), z t splits into and 
z u . The change from J < J c (h) to J > J c (h) is discontinuous. This difference suggests that 
the phase transition is continuous for h = 0 and discontinuous for h > 0. 
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III. ASYMPTOTIC BEHAVIOR OF C(t) 


In this section, we derive the asymptotic form of the correlation function C(t) using the 
previous results for J ^ J c {h). Based on them, we assume the functional form of C(t) for 
J = J c {h). C(t) is defined as the covariance between X(l) and X(t + 1) divided by the 
variance of X (1),'Var (X (1)): 

C(t) = Cov(A(l), X{t + l))/Var(A(l)). 

Through normalization, C(0) = 1. C(t) can be expressed as the difference between two 
conditional probabilities. 

C{t) = Pr (X(t + 1) = 1|X(1) = 1) - Pr (X(t + 1) = 1|X(1) = 0). (2) 

In general, C(t ) is positive for J > 0. 


A. C{t) for J ± J c (h) 


The asymptotic behavior of C(t) depends on (J, h). If J < J c (h), there is one stable fixed 
point at z + and z(t) converges to z + through the power-law relation E (z(t) — z + ) oc f i+_1 126 ]. 


Here, the expectation value E(A) of a certain quantity A is dehned as the ensemble average 
over the paths of the stochastic process. If J > J c (h), there is another stable fixed point at 
Z— Both z± are stable, and z(t) converges to one of the fixed points. The convergence of 
z{t) to z± also exhibits a power-law behavior E (z(t) — z±) oc t l± ~ l 27]. We assume that the 
probability that z(t) converges to one of the z± depends on X(l) and we denote this as 


p±(x) = Pr(^(t) —> z±\X(l) = x). 


For J < J c (h), z(t ) always converges to z + irrespective of the value of X(l) = x, and 
p+(x) = 1 holds. In this case, we set p~(x) = 0. Regarding the asymptotic behavior of the 
convergence of z(t) —> z±, which also depends on X(l), we assume 

E (z(t) -A z±\X(l) = x) ~ W±(x)t l± - 1 
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We write the dependence of the coefficients W±(x) on the value of A"(l) explicitly. Using 
these behaviors and notations, we estimate the asymptotic behavior of C(t) as 


C{t) = Pr (X{t + 1) = 1|X(1) = 1) - Pi(X{t + 1) = 1\X(1) = 0) 

= E(<jl>(«))|A'(l) = 1) - E(«(*(«))|A(1) = 0) 

1 

= ^(—l ) x_1 {E(q(z(t))\x)PY(z(t) —> z+\x) + E(g(A(f))|a;)Pr(A(t) —» Z-\x)} 

X=0 


~ ^(-l) x 1 {(g+ + l+E(z - z+\x))p+{x) + {q- + l-E(z - z_\x))p_(x)} 

X=0 


= 5Z(-1) X 1 {(<?+ + l+W + (x)t l+ l )p+(x) + (g_ + l_W-(x)t l - 1 )p- (a?)} 

x=0 


= ^ [ q+P+(x ) + q_p_(x) + ( l + W + (x)p + (x)t l+ 1 + l_W_(x)p_(x)t l - x )] (-1)* \ 

x=0 


(3) 


Here, we expand q(z) as 

q(z) = q(z± + l±-(z- z±)) ~ q± + l±-(z- z±). 

Given that p+(x) + p~{x) — 1 for x — 0,1, the limit value c = lim^oo C(t) is estimated to 
be 

c = (q+ - ?-)(p+(i) -p+(o))- (4) 

For J < J c (h), p+(x) = 1 and c = 0. As is stable for J > J c (h), the probability for the 
convergence of z(t) to is positive. It is natural to assume that p + ( 1) > p+(0) and c > 0 
for J > J c {h ). 


The asymptotic behavior of C(t) is governed by the term with the largest value among 
{7 + , /_} for J > J c (h). We define l max as 


l 


max 


{ l+ ,, </ < </ c (/i), 

Max{/ + ,/_} , J > J c (h). 


(5) 



We summarize the asymptotic behavior of C(t) as 


C{t) czc + d -t l 1 and l = l max . (6) 

Here, we write the coefficient of the term proportional to t l ~ l as d. If J > J c [h ), the constant 
term c is the leading term. If J < J c {h), the power law term d ■ t l ~ l is the leading term. 
There also exists a sub-leading term to d ■ t /_1 that we do not write explicitly. One reason 
for this is that we do not understand the asymptotic behavior. The second reason is that 
our interest is focused on the value of l. 

B. J = J c {h) 

On the boundary J = J c (h), there are two stable points z + and z t for h > 0. As z t is 
stable, the probability for the convergence of z(t) to z t is positive. It is natural to assume 
that p + ( 1) > p+(0) and c > 0. If h — 0, there is only one stable point at z + = z t = 1/2 and 
c = 0. As l max — It — 1, we anticipate that | C(t) — c\ becomes a decreasing function of In t. 
One possibility is a power-law behavior of In t such as 

C\t) ~c + c'- (lnt)- Q \ (7) 

In the case of the digital model, C{t) oc t~ a with a = 1/2 for p = p c . We denote the power 
law exponent for Inf as a'. 

We derive a' by a simple heuristic argument. At first, we consider the case of h — 0. 
There is only one stable touchpoint at z t , and z{t) converges to it. Eq. (HJ1) suggests that the 
asymptotic behavior of C(t) is governed by E (z — z t \x) as z t is the only stable state. As 
q(zt) = z t and q'{z t ) = 1 at z t , q(z) can be approximated in the vicinity of z t as 

q(z) = -S(z - z t ) 3 + z. 

Here 5 is a positive constant, as z t is stable (FigJSK). The time evolution of E(z — z t \x) is 
given as 

E (z(t + 1) - z t \x) - E (z(t) - zt\x) = ^-yE(g(^(f)) - z(t)\x) ~ “ z t f\x). 
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Here the denominator t + 1 in the middle of the equation comes from the fact that there 
occurs a —h- change in E(z(t)\x) for X{t + 1) G {0,1}. We also assume E((^(f) — z t ) 3 \x ) ~ 
E (z(t) — zt\x ) 3 and the equation can be written as 

^E (z(t) - z t \x) = -jE (z(t) - z t \xf. 

The solution to this shows the next asymptotic behavior 

E(z — z t \x) oc (In 


and we obtain a' — 1/2. 

Likewise, for h > 0, there are two stable states q + and q t . The subleading term in C{t ) 
is governed by E (z — z t \x). We can approximate q(z) in the vicinity of z t to be 

q{z) = S(z - z t ) 2 + z. 

Here 5 is a positive constant (Figj3jo) . If z(t) > z t , z(t) moves toward the right-hand 
direction, on average, and converges to z + . We only need to consider the case z{t) < z t and 
z(t) converges to z t . In the case, E(^(t) — z t \x) obeys the next differential equation. 

^■E (z(t) - z t \x) = ^E (z(t) - z t \x) 2 . 

The solution shows the next asymptotic behavior 

E(z(f) — z t \x) oc (lnhp 1 , 


and we obtain a' — 1. 


C. Numerical check of C(t) ~ c + d ■ t l 1 

We perform the numerical integration of the master equation of the system and check 
the asymptotic forms for C(t) numerically. We denote the joint probability function for 
Y!' s =i X{s) and A"(l) as P(t,n,x i) = Pr(^* =1 A(s) = n,X( 1) = xi). For t — 1, P(l, 1,1) = 


10 


<7(1/2) and P(1,0,0) = 1 — g(l/2). The other components are equal to zero. The master 
equation for P(t,n,x i) is 

Pit + l,n,xi) = q{(n - 1 )/t) ■ P(t,n - l,xi) + (1 - q(n/t )) • P{t,n, xi). (8) 

We impose the boundary conditions P(t, n,x i) = 0 for n < 0 or n > t. Using P{t,n,x i) for 
the case when t <T, we estimate C{t ) for t < T. 



FIG. 4. Plot of c and l in Eqs.© and (fTOl) for h E {0.0,0.5}. □(■): c(l) for h = 0, and o(«): c(l) 
for h = 0.5. The solid (broken) line shows l m ax for h = 0.0(0.5). 


We adopt h E {0.0, 0.5} and J E [0,3], and estimate C(t ) for t < T = 4 • 10 5 . The 
threshold values J c {h) are J c (0) = 1 and J c (0.5) ~ 1.953. We assume the asymptotic 
behavior C(t ) ~ c + c' ■ t for J > J c (0), h = 0 and J > J c (h), h > 0. Using the three 
values of C(t ) at to — T,t\ — T/s, and t 2 = T/s 2 with s = 2, we solve for c, c' and / in the 
following manner. 


I 


C(t 2 ) - C(U) 
^{t^-citoY 
C(ti) — C{t 0 ) 

, 1-1 _ , 1-1 ’ 

L i b o 


c = C(to ) — a ■ t l 0 1 . 


(9) 


For J < J c (0),/i = 0 and J < J c (h),h > 0, we assume the asymptotic behavior C(t) ~ 
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d ■ t l 1 . In this case, we estimate /, c' to be 

^ — 1 — ln s C(ti)/C(t 0 ), 

d = C(t 0 ) ■ tl~ l . (10) 

c and l are plotted using symbols in Fig. HJ We also plotted l max of Eq. (l5]l with solid and 
broken lines. The estimation of / obtained using Eqs. (l9|) and (ITOT) is nearly consistent with 
lmax ■ However, for 1.2 < J < 1.5, h = 0.0 and 2.1 < J < 2.4, h = 0.5, large discrepancies are 
observed in the estimation of l. A possible explanation for these discrepancies is that the 
system size T is not sufficiently large, and therefore the assumption of the aforementioned 
asymptotic behavior does not hold. In particular, as t increases, the sign of c' changes at 
some particular t value, and this negatively affects the estimation of l. We believe that 
these discrepancies can be removed by increasing the system size T. For J = J c (h), the 
aforementioned fitting procedure does not provide accurate results. In the next section, we 
apply the FSS method and clarify the asymptotic behavior of C(t) for J = J c (h). 


IV. FINITE-SIZE SCALING ANALYSIS 


The asymptotic behavior of a system is governed by a temporal length scale, which is 
referred to as the correlation l eng th £. For the definition of £, we adopt the second moment 
correlation length of C(t ) 2l|, 28L which has been adopted in the study of the equilibrium 


phase transition of spin models 


29] and in the percolation theory 30]. We denote the n —th 


moment of C(s) for the period s < t as M n (t ) = X]s=o C(s)s n . The variable t in M n (t ) is 
considered as the time horizon or system size of the stochastic process. The second moment 
correlation length £(t) is defined as £(£) = \JM 2 (t)/M 0 (t). The integrated correlation time, 
also referred to as the relaxation time, r(i), is defined as r(£) = Mo(t). £(£) and r(£) have 
the same dimensions; the length scale in the critical behavior of the system is given by £. 
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A. FSS for J / J c (h) 


We estimate M n (t ) using the asymptotic behavior of C{t) in Eq.( 


t 


pt j-n +1 fn+l 

M n (t) = y ]C(s)s n ~ / C(s)s n ds = c—— + c'- 

J o 


s=0 


n + 1 n + l 


£ t (f) = £(f)/t = y/M 2 (f)/ M 0 (t)t 2 is then given as 


6(f) - 


c = 0 


i (1 + 

3 V T 2cW+2 


, c> 0. 


( 11 ) 


converges to -y/l/3 and + 2) for c > 0 and c = 0, respectively. 

In the case c > 0, we denote the deviation of from the limit vale \J 1/3 normalized by 
the limit value as A£ t . 


A6( () . 


£ t (t) behaves as 


6(f) = 1/3(1 + A&(f)). 


We also describe Tt(f) = r{t)/t as 


T ,(<) ~ C +J l - 


( 12 ) 


We express r 1 using A£ t (t) for c > 0 as 


r,(t) = c 1 1 + 


(13) 


converges to c as A£ t (f) converges to zero. 

We define the scale transformation of the system as the change in the time horizon of 
the system t —> st with the scale factor s. We denote the scaling function for a long-term 
observable A{t) as /^, which is defined as 


A(st ) 


— /a (6(f)) + A/a(f). 
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Here, we include the correction to the scaling term as AA/^(f) is of the order t 


and oo is a correction-to-scaling exponent 


m 


When C(t) exhibits power-law behavior, the system is scale invariant 21]. Furthermore, 
linit-Kjo £ t (£) is constant and = 1. Afor c > 0 is given by 

A/ft= w ~ 1 ~ (S,_I “ 1)A? " 


As A£t oc t l 1 , oo = l — 1 for 


If c > 0, f Tt = 1 as lim^oo r t = c, then 


Mn = 


Tt(st) l + 2 , i_i 


nit) 


l ~ 


/-l 


1)A6- 


and oo — l — 1. If c = 0, r t oc r 1 and 


In 


; Jr t = lim In, —— 

t—loo 


-2 


= l - 1 = 


3-6 

er 2 -i' 


(14) 


As l < 1, under the scale transformation t —> st, Tt(st) = s l 1 Tt(t) oc s l In the limit 
s — y oo, Tt{st) decreases to zero. 


B. FSS for J = J c (h) 

We consider and verify Eq.Q by studying the finite-size scaling relation of the system. 
We estimate M n (t ) to t n+1 (lnt) _Q,_1 as 

/ i fn+1 rt 

C{s)s n ds ~ c-f d / s n (ln s)~ a 'ds 

n + 1 7 


^n+1 ^n+1 


~ C 


+ c-- (Inf) “ 1-|--(Inf) 


n + 1 n + 1 


a 


v-l 


n + 1 
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Here, we use the next expansion in powers of 1/Inf, which is obtained by partial integration. 


s n (lns) a 'ds = 


1 r +1 (lnf)-“'+ “ 


n + 1 


ra+1 


71+1 


t n+ 1 (lnt)-° + 


a 


77+1 


f n+1 (lnf) _a ' _1 + 


f n+1 (lnf)- Q ' + 


(77 + 1 )' 


a 


(77 + 1 )' 


s n (lns 

a'(a' + 1) 

(n + l ) 2 


, —Q—1 


ds 


s n {lns 


a '—2 


ds 


r+^lnf)-"'- 1 + 0(f n+1 (lnf)-"'- 2 ) 


n+1/ 


In the following, we omit the parentheses (,) in (Inf) x when we write the power of Inf. 
£t{t) is 

|v^( 1_ ¥ lnrl ) ’ c = 0 

lv^( 1 _ ^( lnrQ ' _ 1 “i lnt ' 2a ')) > c >°- 

£t converges to the same value ft 1/3 in both cases, c = 0 and c > 0. However, the 
convergence speeds differ. The coefficients and exponents of the power of Inf of the sub¬ 
leading terms are given as B and f3 ', respectively. 


6 - V g(l - B Inf P) = 


(i+A m) 


(15) 


Using the second equality, we define A£ t . f3' = 1 for c = 0 and ft = Min(l + a', 2a') for 
c > 0. 


r t (t ) is estimated to be 


r t ~ c + d In f . 


(16) 


As converges to ftl/3, the system is scale-invariant and = 1. For c = 0, we can 
use the scaling relation for r t in order to estimate a’. Because T t is a function of Inf, we 
consider the exponential scaling transformation f — y t s with the scale factor s. 


Bit 8 ) = , 

nit) 


(17) 


For c > 0, as f Tt = 1, it is necessary to study the finite-size scaling correction for r t . 

- 1 ~ -(s" Q ' - 1) In t~ a '. (18) 

nit) C K J y J 
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We can estimate cx' and c'/c using Eq. (lT8j) . In addition, we can estimate (3 ' using £ f . 


6(0 

6(0 



i)A6- 


( 19 ) 


If cx' = 1, a convenient extrapolation formula for c is available. As = — -^(lnf)^ 2 , we 
can express c'/c using and Ini. We obtain the extrapolation formula for c by solving 
Eq. ffTfij) as 


n(t) 

1 — 4A £ t (t) Ini 


( 20 ) 


V. NUMERICAL STUDIES OF FSS AND C(t) FOR J = J c (h) 

We check the assumption C(t) ~ c + c'lnt~ a ' by studying fini te-size scaling. First, we 
estimate cx' for h = 0 using Eq. ffTTj) and the data for t s = 3 x 10 6 and s G [1.0,1,1]. 
Figure [5](a) shows the plotted results. It is evident that the data lies on the curve s -1 / 2 and 
cx' = 1/2. 

For h > 0, we apply Eq. (fl8]) in order to estimate ex'. We estimate T t {t s )/T t {t ) — 1 for 
s = 1.01,1.001 and divide it by s~ a ' — 1 using the data for h G {0.1, 0.3, 0.5} and t s < 5 x 10 6 . 
We choose ex' so that two data sets with different s lie on the same curve as a function of 
1/lnf. Figure E][b) shows the double lnarithmic plot of the obtained results. We choose 
ex' = 1 and two data sets for each h that lies on the same curve. If the curve obeys Eq. (lT8lh 
the slope should be equal to 1 and then, we can estimate c'/c by fitting with c'/c In t. 
However, even for t s = 5 x 10 6 , the slope of the curve is not equal to one for h = 0.1. We 
suppose that the system size t s = 5 x 10 6 is not large enough for Eq. IjTSj) to hold in this case. 

In order to verify that cx' = 1 for h > 0, we use Eq. (fT9l) . As /3 1 = Min(l + ex', 2cx'), j3' = 2 
for cx' = 1. We use the estimate of £t(t s )/^(t) — 1 with the same data set and s = 1.01,1.001. 
We divide it by s~ 2 — 1. Figure [5](c) shows the plotted results. The two data sets for each 
h lie on the same curve; this agrees with our estimation /3' = 2. Furthermore, as t becomes 
large, the data set converges to A£. These results are consistent with Eq. dlOl) and confirm 
that cx' = 1. 

We estimate the fitting parameters c, c' using the numerical results of C(t ) for t < 5 x 10 6 . 
After c is subtracted from C(t), C(t ) — c behaves as c'(lnf) -1 for c > 0. Figure [6] shows a 
double lnarithmic plot of C(t)—c vs. In t for J = J c {h). It is difficult to verify the power law 
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dependence of C(t) — c on Inf using the narrow range of Inf. The data lies on the straight 

lines. 



FIG. 6. Plots of C(t) — c vs. Inf for J = J c (h ) and h € {0.0, 0.1,0.3, 0.5}. t < T = 5 x 10 6 (3 x 10 6 ) 
for h > 0 {h = 0). For h = 0.0, we plot the fitted results using C(t) = c'lnf -1 / 2 with a fitting 
parameter of d. For h > 0.0, we plot the fitted results using C(t) = c + c'lnt -1 with two fitting 
parameters c, c!. 


A. Extrapolation of c 

We now estimate the limit value c = lim^oo C(t) for J = J c {h). A simple method is to 
use the fitted result with the assumption that C(t) = c + c' Inf -1 . Another method is to use 
Eq. (l20j) . We summarize the results of these methods in Table |H 

TABLE I. The columns in the table denote the following: an estimation of c for h € {0.1,0.3,0.5} 
and J = J c (h). h in the first column; C(T — 1) in the second column; the estimate by fitting 
in the third column; tj(T) in the forth column; the estimate using Eg. (12011 in the fifth column. 

T = 5 x 10 6 . _ 

h C(T — 1) c (Fitted) r t (T) c (Extrapolated) 


0.1 

0.558 

0.504(1) 

0.560 

0.514 

0.3 

0.638 

0.582(2) 

0.640 

0.593 

0.5 

0.659 

0.604(2) 

0.661 

0.615 


Both C(T — 1) and T f (T) provide slightly larger estimates of c than the fitting and 
extrapolation methods. The fitted and extrapolated values are approximately 10% smaller 
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than the estimates using C(T — 1) and r*(T) for T = 5 x 10 6 . An extremely slow convergence 
of C(t) to c is observed. 


VI. UNIVERSALITY CLASS 

The system shows a continuous phase transition for h = 0. As we have seen, the analysis 
does not depend on the precise form of q(z). As far as q(z) has ^-symmetry, g( 1 — z) — 
1 — q(z) and y = q(z ) are tangential to y — z at z t — 1/2, we can assume a similar continuous 
phase transition. In order to discuss the universality class of the phase transition, we compare 
the scaling properties of C(t) for two models. One model is the logistic model which adopts 
the q(z) in Eq.fll]). For h = 0, J c (0) = 1. Another model adopts the next q r (z) with three 
parameters r, p e [0, l],g* 6 (1/2,1][? ]. 

q r (z) = (1 - p) ■ Q* + P ■ n r (z) 

r 

■K r {z)= Y, rC s -Z S {l-Zy~ S (21) 

s=(r+l)/2 

Here r C s is the binomial coefficient and r is a odd number greater than three, r G {3, 5, 7 • • • }. 
In the limit r —» oo, q r (z) reduces to that of the digital model, (1 — p) ■ g* + p ■ 9{z — 1/2). 
This model corresponds to the mean-field approximation of the model, where X(t) chooses 
the majority of r randomly chosen previous variables with a probability of p. For g* = 1/2, 
q r (z) has ^-symmetry and the threshold value p c {r) is determined by the condition 1 = 
q' r (z t = 1/2) = p c {r) ■ 7t'(1/2). p c is explicitly given as 

Pc{r) =- - -. 

p c {r) is a decreasing function of r and lim r _ ) . 0O p c (r) = 0, which is compatible with p c = 
1 — l/2g* = 0 of the digital model with g* = 1/2. 

At J = J c { 0), C(t) obeys a power-law of In t as C(t) ~ c / lnt^ a, with a' = 1/2 for the 
former model. Below J c (0) = 1, C(t) oc t 1 ^ 1 with l — l + — J. We set A J = J c — J = 1 — J. 
The expression for the exponent of C(t) oc 1 1 ~ l is given by 

1 - l = AJ , J < J c ( 0) = 1. 
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For J > 1, C{t) — c oc r 1 with / = /+ = /_ = q'(z + ). As in the case of the estimation 


24], we solve £+ = q(z + ) with the 


of the critical exponents for the mean-field Ising model 
assumption A J = J — 1 << 1. We obtain z + — 1 = \/3A J and estimate l as 


1 - l = -A J , J > J c (0) = 1. 


As J approaches J c (0) from below and above of J c , l — l approaches 0. At J = J c , C\t ) 
obeys a power-law of In t. As = e~d _i ' lnt , we can regard 1/(1 — /) as the ’’correlation 

length” while assuming C{t ) as function of In t. We assume the phenomenolnical scaling 
ansatz for C{t) to be 

C(t) = lnt - " g((l — /) lnt). (22) 


In the ansatz, In t dependence of C{t ) is scaled with 1/(1 — /). g(x) is the universal function 
and is finite at x — 0. For J = J c and l = 1, C(t) = g(0)lnf~ Q/ with g(0) = c!. In the 
limit x —> oo, in order to compensate the In tr a ' term, g(x) should behave as g(x) oc x a> for 
J > J c ( 0). Then C{t) behaves as lim^ 00 C(t) oc (1 — l) a ' = A J“\ The critical exponent 
for c oc A J 13 coincides with a'. The exponent v\\ is defined for the divergence of 1/(1 — l) as 


1/(1-/) oc AJ^ll. 


As (1 — /) oc A J, we obtain v\\ = 1. The scaling relation /3 — a 1 ■ v\\ holds. 

For the latter model with q r (z), l = p ■ q' r (z t = 1/2) for p < p c (r). As 1 = p c ■ q' r ( 1/2) 
holds, the correlation length 1/(1 — /) is estimated to be 


<(1/2 )(p c (r)-p) 

and diverges as 1/(1 — /) oc A p~ u ^ with v\\ = 1. Here, we define A p = p c (r) — p. For 
p > Pc, we assume A p = p — p c (r ) << 1 and estimate z + by solving z + = q r (z + ) up to 
0(Az 3 ) = 0((z + — 1/2) 3 ). One may then show that 1 — l = q' r {z + ) oc A p, and we obtain 
z/|l = 1. If we assume the scaling form for C(t) to be C(t) ~ ln^ _a, • g(lnt ■ (1 — /)) and define 
/3 as c oc A pP, we obtain /3 — a' ■ v\\. 

If the two models share the same value for a this suggests that they are in the same 
universality class. 
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A. Numerical calculation of g(x) 


We estimate the universal function g{x) assumed in Eg. (1221) numerically. For the former 
logistic model with h = 0, we estimate that a' = 1/2. J c (0) = 1 and we estimate C(t) 
for t < 4 x 10 5 and 2/3 < J < 3/2. For the latter model with r — 3 and g* = 1/2, 
p c { 3) = 2/3. We estimate C(t) for t < 4 x 10 5 and 4/9 < p < 1. Using data for C(t) 
between 10 4 < t < 4 x 10 5 , we determine g{x ) to be 

g(x) = lnt 1 / 2 • U(t) , a; = (1 — l) lnt. 

g(x) should be smooth near x = 0 and g(0) = d For a sufficiently large J, c ~ 1 and l = 0. 
g(x) should behave as x 1 ^ 2 for sufficiently large x values. For J < 1, g{x) should decrease 
exponentially. 



FIG. 7. Plot of lnt 1 / 2 • C(t) vs. (I — l)lnt. 

Figure [7] shows the results of this analysis. The thick continuous and thick dashed lines 
indicate the results for J > J c ( 0) and p > p c ( 3), respectively. As can be clearly seen in 
this figure, the data obtained for different J and t values and for different p and t values lie 
on two curves, which represent g(x) in the phase with c > 0 for both models. For large x, 
g{x) — x 1 / 2 . The thin continuous and thin dashed lines indicate the results for J < J c (0) 
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and p < p c ( 3), respectively. The data lie on two curves, which represent g{x) in the phase 
with c = 0 for both models. g(x) can be seen to decay exponentially. The results indicate 
that the scaling ansatz in Eq . (122 1) holds with a' = 1/2. 


VII. SUMMARY AND NOTES 


In this study, we analyzed the asymptotic behavior of the normalized correlation function 
C(t) for a generalized Polya urn. The probability q(z) of adding a red ball to a specific 
proportion of red balls z is q(z) = (tanh[J(2z — 1) + h\ + 1)/2. There are three domains in 
( J,h ): J < J c (h ), J > J c {h), and J = J c {h), where for J = J c {h), q(z ) becomes tangential 
to z. The limit value c = lim^ 00 C'(i) is the order parameter for the phase transition. If 
J < J c (h ) (> J c (h)),c = 0 (c > 0). C(t) ~ c + a ■ t 7-1 with c = 0 (> 0) for J < J c (J > J c ), 
and l is the (larger) value of the slope(s) of q(z) at stable fixed point(s). Through FSS 
analysis, we evaluated the asymptotic behavior of C(t) for J = J c [h ) as C(t) ~ c+c' ln(t) _a \ 
For h = 0, c = 0 and a' = 0.5. For h ^ 0, c > 0 and a' = 1. The system shows a continuous 
phase transition for h — 0. C(t) behaves as C(t) = ln/V Q, g((l — /)lnt) with a universal 
function and we numerically estimated g(x). The scaling relation j3 = ol'v m holds true among 
the critical exponents with a' = 1/2, mi = 1 and /3 = 1/2. We also studied the scaling of 
C(t) for q r (z) given in Eq. d^ID with z * = 1/2 and r = 3. We showed that the scaling ansatz 
C(t) = logt~ a 'g((l — l ) logt) holds for a' = 1/2. 

We note several key points and future challenges. The first of these is related to the 
relationship between the non-equilibrium phase transition studied in this work and equilib¬ 
rium phase transition of the mean-field Ising (MFI) model for h = 0. As is well known, 
/3 = 1/2 for the latter model and both types of phase transitions give the same values for 
(3. A crucial difference lies in the behavior of C(t). In the equilibrium phase transition, the 
memory of the value of X(l) should disappear for J < J c (0). However, in the case of the 
non-equilibrium transition, C(t) decays according to a power-law of t and continues to exist 
for finite t. Furthermore, for h > 0, phase transition does not occur in the equilibrium case 
and z(t) converges to z + . In the non-equilibrium case, the probability for z(t) to converge 
to is positive for J > J c (h). In our previous study, we controlled the length of the mem¬ 
ory r, and also, X(t + 1) depended only on recent r variables X(t — r + 1), • • • ,X(t) [25]. 
The Polya urn process corresponds to the case where r = t. We have shown that with the 
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logarithmic increase of r oc Inf, the non-equilibrium phase transition for h > 0 disappears 
and z{t) always converges to z + . We believe that it is possible to understand the relation 
between the non-equilibrium and equilibrium phase transitions by controlling the increase 
of r. We assert that if the increase in r is infinitely slowly and the variables are completely 
equilibrated among the recent r variables, the non-equilibrium phase transition reduces to 
the equilibrium phase transition. 

A problem for the future is related to the derivation of a' and g(x). Our study only 
provides numerical results and a heuristic derivation of cd; A more rigorous mathematical 
treatment appears to be necessary. As the model with q r (z) shares the same value of a' = 1/2 
for r = 3, and the heuristic derivation of a' only uses the approximate form of q{z) in the 
vicinity of the touchpoint, the universality class of the continuous phase transition of a 
generalized Polya, where q(z) has ^-symmetry and q(z) becomes tangential to the diagonal 
at z t — 1/2, is described by a' = 1/2. 

It is also important to verify these results u sing experimental data. Information cascade 


experiments could be one possible solution 
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19]. As the empirically estimated q(z) does 


not exhibit ^-symmetry at p — p c , the results of this paper suggests that the transition 
in the experiment is discontinuous. The system size T is severely limited in laboratory 
experiments. Therefore, web-based experimental systems should be developed in order to 


study the asymptotic behavior of the system 
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Appendix A: g(x) for the ’’digital” model 


We derive the universal function g(x) for the exactly solvable digital model using the 
results given in Ref. [21]. We adopt q(z) = (1 — p) ■ q* + p ■ 0(z — 1/2). The model shows 
a continuous phase transition at p — p c (q*) = 1 — 1/2g* for g* > 1/2. C(t) behaves as 
C(t) — b(q„)t^ 1 ^ 2 at p — p c (q *). The limit value c(g*,p) = lim^oo C(t) is a continuous 
function of g*,p and becomes positive for p > p c (q*)- 


g(x),q*=0.6 



FIG. 8. Plot of C(t)/b{q*)t~ 1 / 2 vs. t/£(q*,p) with empty symbols indicating p > p c (q*) = 1/6 and 
filled symbols indicating p < p c {q*)- We adopt g* = 0.6 and 10 4 < t < 10 5 . The lines show the 
results of Ea. (|A3j) . 


We assume the scaling form for C(t) as 

C{t) = b(q*)t~ 1/2 g(t / £(q*, p)). 
fe(g*) and £(g*,p) are defined as 


b{q*) 

£(<?*, p )" 1 


[8 { 2q.-l \ 

V VT \4 g* -l) 

— In a/4(p + (1 — p)(l — g*))((l — p)q*) 


(Al) 


(A2) 
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g{x) is then given as 


9{x) = 


\J 47 TX + ^ f x °°u 3/2 e U du , p>Pc(q*), 
/“ u- 3/2 e~ u du , p<p c (q *). 


(A3) 


For the derivation of g(x), we use the explicit form of C(t) for g* — 1 given in Ref. 21]. 
We expand £(g*,p) and c(g*,p) around p = p c (q*) and take the limit p —* p c (q *) by fix¬ 
ing x = t/£(q*,p). For the general g* > 1/2, we estimate b(q*) using the expressions for 
c (q*?p)^(q*iP) an d the assumption in Eq. flAlIh For sufficiently large x, g(x) — yJ^KX. As 
lini/^oo C(t) = c(q,p ), for p ~ p c (q*), we obtain 


c(q*,p) = b{q*)/y/€{q*,p). 


We can estimate b(q*) by 


b(q*) 


V £(q*,p) c (q*,p) 

Inn —- -= - 

p->Pc(q*) V47T 


We have derived g{x) exactly only for g* = 1. We are able to check g(x) for g* < 1. We 
can estimate C(t) for t < 10 5 and g* = 0.6 by numerically integrating the master equation 
for the system. Figure [8] shows the results of this process. The empty symbols indicate the 
results for p > p c (q*)- We have adopted the value of p in the vicinity of p c (q*) = 1/6 and 
te[ io 4 ,4 x 10 5 ]. As can be clearly seen, the data obtained for different p and t values lies 
on the curve of Eq. (1A3[) . The Filed symbols indicate the results for p < p c {q*)- The curve 
of Eq. (lA3l) describe the data. 
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